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ABSTRACT 

Structure of a quasi-stationary stellar cluster is modelled assuming that it is embed- 
ded in the gravitational field of a super-massive black hole. Gradual orbital decay 
of stellar trajectories is caused by the dissipative interaction with an accretion disc. 
Gravitational field of the disc is constructed and its effect on the cluster structure 
is taken into account as an axially symmetric perturbation. Attention is focused on 
a circumnuclear region (r < 10"* gravitational radii) where the effects of the central 
black hole and the disc dominate over the infiuence of an outer galaxy. It is shown 
how the stellar system becomes gradually flattened towards the disc plane. For certain 
combinations of the model parameters, a toroidal structure is formed by a fraction of 
stars. Growing anisotropy of stellar velocities as well as their segregation occur. The 
mass function of the inner cluster is modified and it progressively departs from the 
asymptotic form assumed in the outer cluster. A new stationary distribution can be 
characterized in terms of velocity dispersion of the stellar sample in the central region 
of the modified cluster. 
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1 INTRODUCTION 

We discuss the impact of gaseous environment on the long- 
term stellar motion in central regions of active galactic nuclei 
(AGN). We adopt the view of galactic nuclei as composite 
systems with stars surrounding a super-massive black hole 
(with a typical mass of W^Mq <M.< Mq) and forming 
a rather dense cluster (Begelman & Rees 1978; Norman & 
Silk 1983). Central number density of the cluster stands as 
one of the model parameters that can reach values of the 
order of n* ~ 10^ -ir lO^pc"^. We consider the stars of the 
cluster as satellites orbiting the central black hole. We also 
take into account the influence that the diluted interstellar 
environment exhibits on the motion of stars on time-scales 
exceeding 10* dynamical periods. This allows us to address 
an important question whether or not the presence of such 
medium could increase the concentration of stars in certain 
regions of the nucleus. 

Although we consider the region of direct gravitational 
influence, where the gravity of the central black hole dom- 
inates over non-gravitational forces, one should not neglect 
the impact of perturbing effects on individual stellar trajec- 
tories and the overall cluster structure. We assume that it is 



the gas in an accretion disc which acts as the main pertur- 
bation. The disc medium affects the motion of stars either 
directly by collisions between stars and gas particles at the 
point of transition through the disc slab (the case of inclined 
trajectories), or vicariously, via coupling between the stars 
and the disc medium (embedded trajectories). These pro- 
cesses will be discussed in combination with the effects of 
disc gravity. 

Our present contribution follows a previous paper 
(Vokrouhlicky & Karas 1998) where the long-term evolu- 
tion of stellar orbits was studied including the case of an ax- 
ially symmetric perturbation imposed on the gravitational 
fleld to mimic the presence of a disc. It was shown that the 
orbits experience abrupt changes of eccentricity and incli- 
nation even if the disc mass is a tiny fraction of the central 
mass, Md < 10~"'M,. In that paper properties of individual 
orbits were investigated. A simplified approach to the evolu- 
tion of the central cluster was adopted in Ranch (1999) and 
Karas, Subr & Slechta (2002). Stellar encounters were ex- 
plored in this context by Vilkoviskij & Czerny (2002), using 
a semi-analytical approach to the dense and flattened clus- 
ter. Here we add an astrophysically more realistic descrip- 
tion of the self-gravitating disc, especially of its gravitational 
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field that superimposes with the central field, and we study 
dynamics of the cluster in terms of its overall characteristics 
(the shape and velocity dispersion). 

The main new ingredients of the present investigation 
are threefold. Firstly, a more elaborate description of the 
disc accounts for different processes and corresponding time- 
scales that dominate the system evolution at different dis- 
tances from the core. Orbit behaviour is governed mainly 
by the overall radial profile of the disc (especially its surface 
density and vertical thickness) and by stellar characteristics 
(mass and size). The interplay between these principal pa- 
rameters determines the hydrodynamical mechanism driv- 
ing the evolution - collisions with the disc, density waves, 
or formation of gaps in the disc. 

Secondly, we introduce the distribution of stellar masses 
in the cluster, A/'(M«). The cluster is supplied with a 
Salpeter type initial mass function, but the resulting mass 
function is eventually different because of selective orbital 
decay. For our purposes, stars are characterized phenomeno- 
logically by their effective column density Et=Af,/ (vri?*)^ 
as a key paramater that rules this effect. 

Thirdly, we consider the distribution of stars in phase 
space. Simulations start with a spherically symmetric and 
stationary cluster whose dynamics is dominated by the cen- 
tral black hole mass (BahcaU & Wolf 1976, BW hereafter). 
This solution describes a gravitationally relaxed system of 
stars near a black hole with phase-space probability den- 
sity /0(a) oc a^''* (/o(a) da is the number of stars in a shell 
with semi-major axis a and width da). We assume that the 
modified distribution approaches /o(a) far away from the 
core, but it can be modified in the inner parts of the system 
where also the importance of the disc increases. Other self- 
consistent models have been discussed for a cluster with a 
central black hole (e.g. Young 1980; Quinlan, Hernquist & 
Sigurdsson 1995; Takahashi & Lee 2000; Alexander 2003) for 
which the radial density profile turns out to be different but 
also a power-law type at large radii. Furthermore, several 
authors (e.g. Ranch 1995; Freitag & Benz 2002; Alexander 
& Hopman 2003) have developed numerical approaches tak- 
ing stellar collisions into account in a dense cluster near a 
black hole. 

The paper is organized as follows. In the next section 
we describe the ingredients of the model and time-scales in- 
volved. In Section |3 we show a modified quasi-stationary 
stellar distribution that emerges in the inner regions of the 
cluster. We discuss transition regions where the stars are 
gathered and we also examine predicted velocity dispersions, 
mean spectral-line profiles and the modified mass function. 
Some of these properties may have direct observational con- 
sequences. Finally, limitations of the model are discussed 
and conclusions are presented in Section 2] 



2 THE MODEL 

2.1 Components of the model 

Our exploratory model is based on a presumably obvious 
assumption of dissipative processes prevailing in the center, 
where the disc surface density is high, while gravitational 
relaxation is supposed to be in control on the outskirts of 
the system behind the outer edge of the disc. The adopted 



approach contains simplifications that help us to proceed 
in the situation when Af-body or equivalent methods are 
still a challenge, especially if the dissipative environment, 
different types of stars and non-spherical geometry of the 
system are to be taken into account. The model consists of 
three components: 

(a) A super-massive black hole. Its infiuence is approx- 
imated by Newtonian gravitational field of a point mass 
Af, . Relativistic effects are introduced only by adopting 
Rs=GM./c^ ~ 1.5 X 10^^ Mgcm {Ms=M./10^ Mq) as the 
smallest natural length-scale that appears in our consid- 
erations. The corresponding dynamical time-scale is much 
shorter than all other time-scales in the model and it does 
not stand explicitely in equations of the long-term evolution 
hereafter. 

(b) An accretion disc. In order that computations be well 
defined and easier to compare, we consider the standard 
Shakura-Sunyaev disc (denoted as "SS"; e.g. Frank, King 
& Raine 2002) and a model of weakly self-gravitating disc 
(Goldreich & Lynden-Bell 1965; Paczynski 1978; Shlosman 
& Begelman 1989). We assume that approximation of ver- 
tically averaged quantities is sufficient for our analysis and 
we characterize the disc by surface density ^d{r) and geo- 
metrical semi-thickness H{r) (functions of radius in the disc 
mid-plane z = 0). Other variables are also relevant for star- 
disc interactions, e.g. radial velocity of mass transport in 
the disc, temperature of the medium and speed of sound. 
Gravitational field is determined numerically. With the disc 
we associate time-scale tinci, i.e. the interval on which stellar 
orbits are inclined into the plane by the dissipative interac- 
tion. 

(c) A nuclear cluster. It is characterized by the distribu- 
tion function /* that conforms initially to a spherically sym- 
metric and stationary state in the central potential well. In 
our computational scheme, the cluster is divided in two dis- 
tinct subsets that are treated separately. In the inner cluster, 
the stellar system becomes oblate and anisotropic due to the 
dissipative interaction with the disc. Efficiency of this inter- 
action depends on the parameters of individual stars and, 
therefore, it is sensitive to their masses, orbital inclinations, 
as well as to properties of the dissipative environment. Grav- 
itational relaxation is neglected in the inner cluster, assum- 
ing that ti ^ tinci- The outer cluster, on the other hand, 
does not evolve in our computations and its structure is 
determined by gravitational relaxation as a dominating pro- 
cess. The outer cluster acts as a reservoir from which fresh 
stars are injected inwards, maintaining the overall steady 
state. 



2.2 An accretion disc with effects of self-gravity 

We adopt three fiavours of the standard model: (i) radia- 
tion pressure dominated standard (SS) disc; (ii) gas-pressure 
dominated SS disc with electron scattering; and (iii) the 
same as the previous case but with free-free opacity (Frank 
et al. 2002). Radial dependence of the disc variables is ex- 
pressed by power-law in all those cases, in particular, the sur- 
face density can be written in the form Ed/E* oc [r/RgY . 
The values of index s are given in Table Q for different cases 
(see also Karas & Subr 2001) and we will introduce also 
non-power law profiles later in this section. 
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Figure 1. Gravitational field of different models of a disc is shown in comparison with the central field of the black hole. Ratio 
of the corresponding field components was computed and the resulting contours plotted (in azimuthal section). The ratio of vertical 
{z) components is shown in the top row (i.e. disc/^z, centre = const), while the ratio of radial (cylindrical R) components is in 
bottom (Ffl.diaeZ-F'fl, eentre — const). Obviously, near the equatorial plane the vertical component of the disc gravity dominates. Three 
representative types of discs have been selected for illustration, as indicated above the frames. Contours are drawn in logarithmic spacing 
and the values refer to = 1. 



Table 1 . Values of the power-law index s characterizing the radial 
dependency of the disc surface density Ti^(r). 

Model (i) (ii) (iii) (iv) (v) 

s parameter 3/2 -3/5 -3/4 -15/7 -5/3 



The standard scheme was supplemented by consider- 
ing self- gravity according to Collin & Hure (1999, 2001) 
and Hure et al. (1994, 2000). Based on Toomre's crite- 
rion (Q ~ 1) for onset of self-gravity, AGN accretion 
discs are affected beyond critical radius ~ 10''J?g (i.e. 
0.01 pc, depending on model parameters). Farther out, at 
> 1 pc from the center, the disc mass can reach a size- 
able fraction of the central mass. For example, the mass 
ratio of the disc (ii) with respect to the central mass is 

H = Ma/M. = O.aMg/^Mo/i^Q^t'''^ (i?d/10''_Rg)^^'^ where 
we J?d is the disc outer radius, Mo.i = M/O.lMEdd, 
Qo.i = a/0.1. Similarly, for case (iii) we have ^ — 

0.1Ml^^M^{^°ao t^'^ {Rd/W''RgY^'^, suggesting that the 
disc attraction is strong enough to influence the motion of 
stars independently of the details of a particular model. Self- 
gravity manifests itself by gathering material to the mid- 
plane, thereby narrowing vertical thickness. 

We considered the case of marginally stable self- 
gravitating discs. Two examples, denoted as (iv) and (v) 
hereafter, are characterized by Q ~ 1 and they differ by 
assumption about elemental abundances: solar metalicity is 
adopted in the case (iv) while it is set to zero in the case (v) . 
Because of differences in opacities, different density profiles 



and values of the total disc mass are obtained. Typically, 
in the outer region of a marginally stable disc, density de- 
creases faster compared to the standard case. For a detailed 
discussion of the adopted models, see Hure (2000) and ref- 
erences cited therein. 

Finally, we considered composite discs, by which we 
mean a combination of simple models. For example, case 
(ii)-l-(iv) is constructed by concatenating two corresponding 
basic models. Notation indicates that the standard SS disc 
(ii) is assumed at small radii and matched to a marginally 
self-gravitating disc (iv) at large distance from the centre 
(density and other variables are continuous at the transition 
radius). This approach reflects an obvious fact that simple 
models can be valid only for a limited range of radius and 
provides setup for the model that is rich enough to capture 
different possibilities. 

Figure illustrates three examples of the structure of 
the gravitational field. To this aim we plot the ratio of field 
components that correspond to the disc and to the cen- 
tral black hole (contribution of the cluster stars is not in- 
cluded in this figure). Different frames reveal the difference 
between the case of Kuzmin's disc (employed in Vokrouh- 
licky & Karas 1998) and two typical examples of our discs, 
computed numerically for this paper. Even if the gravita- 
tional field of the disc is perturbative {fi <^ 1), it turns out 
that the vertical component of the disc field may success- 
fully compete with that of the centre in certain regions near 
the disc plane. This is essential not only for the disc internal 
structure but we can also expect stellar paths to be visibly 
dragged (on the long term), especially those which follow a 
low-inclination orbit, i.e. almost embedded in the disc. 
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The gravitational potential and the field due to disc 
are computed according to Poisson's integrals in which in- 
tegrations are carried over volume occupied by the g 
suming axial symmetry. We used the splitting method for 
accurate evaluation (Hure 2000). Furthermore, in order to 
reduce duration of runs, we implemented a dynamical grid 
in the Poisson solver and we pre-generated the field compo- 
nents in a Cartesian grid covering a rectangular region of 
1 pc X 1 pc. In some regions, close to the disc, the field lines 
are very distorted. We found it convenient to define addi- 
tional sub-grids which provide fine resolution where needed. 
Field values between mesh points were obtained by bilinear 
interpolation. Hereafter, we will illustrate results for discs 
(iii) and (ii)-l-(iv), i.e. the same as in Fig.0 These two cases 
appear sufficient in order to demonstrate various stages of 
the modified cluster. The models share the same black hole 
mass. A'/, — 10*Mq, as well as other common parameters 
of the disc: = IQ'^R^, M = CIA/e, and a = 0.1. Corre- 
sponding mass ratio is /i = 1.2 x 10~^ and ^ = 2.9 x 10~^ 
for discs (iii) and (ii)-f(iv) respectively. 



2.3 Star— disc interactions 

Stars passing through the disc may create gaps and holes in 
some regions, but we assume that the disc is not removed 
completely, implying that the accretion rate cannot be too 
low and the cluster excessively dense. We adopt the approxi- 
mation of impulsive (i.e. instantaneous) interaction between 
crashing stars and the disc (Syer, Clarke & Rees 1991; Karas 
& Vokrouhlicky 1994). In this way, stellar trajectories slowly 
evolve and several phases can be distinguished in the course 
of orbital evolution. 

During the initial phase, as long as the stars are crashing 
on the disc with non-zero inclination, the following combi- 
nations of orbital elements a, x = cosi and y = (1 — e^) 
remain constant: 



Ci = ay (I + x) and C2 



{l-y)il+x)' 
1 - X 



(1) 



(see Subr & Karas 1999). These two conserved quantities 
can be employed to determine the orbital evolution along 
with a simple estimate of tinci, 



tincl — toMs 



ao 



S0 \R, 



3/2- 



(2) 



Here, various model-dependent factors were absorbed in the 
term to=to{ao,eo,io) that varies only slightly with ao. For 
example, in cases (ii), . . . , (v), and for eo = 0, io = 7r/2, we 
find to = 0.92, 0.22, 1.89 x 10"'^ and 9.40 x 10"", respec- 
tively (these numbers were determined numerically).^ Be- 
cause there is considerable uncertainty in the actual value 



^ It can be shown that Ci and C2 are strictly conserved for or- 
bits with pericenter argument ui = n/2, but their variation for 
arbitrary uj is only weak. In the case (i), to would come out very 
small compared to other cases because S(j increases with R. How- 
ever, this is the radiation pressure dominated disc, which can be 
relevant only in the innermost regions, and so it should not be 
used to describe the whole process of inclining the trajectory into 
the disc. We avoid such behaviour by using composite models of 
the disc as described above. 



of the disc surface density and since the cluster consists of 
various stellar types, ranging from compact stars to giants, 
it is easy to realize that tinci must span an enormous range 
over seven orders of magnitude in eq. Q. 

Inclination and eccentricity decay exponentially at late 
stages when i oc exp[— t/r], e ~ \/C^/(4i), while the semi- 
major axis slows down and eventually reaches its asymptot- 
ical value flf = Ci/4. When inclination drops to a critical 
value (given by tan i^H{r)/r), the regime of direct collisions 
goes over to the regime of stars entirely embedded in the 
disc. Since this moment, we further distinguish two modes 
of inward migration through the disc medium. If the stellar 
mass is large enough to create a gap (e.g., Lin & Papaloizou 
1986), the star becomes tidally coupled with the disc and it 
sinks to the centre at roughly the same rate as the material 
of the disc. Otherwise, excitation of density waves dominates 
and causes radial migration to be faster, typically by order 
of magnitude (Ward 1986; Artymowicz 1994). 

In Figure 121 we show a particular case of an orbit evolv- 
ing with the disc (iii) and A/* = 3Mq, ensuring that the 
star successively enters all three modes of migration. Starting 
from ao — 10"^ -Rg, eo = 0.45 and xq = 0.54, the semi-major 
axis converges to at = 0.47ao along with the simultaneous 
decay of eccentricity and inclination. The initial phase of 
star-disc collisions is followed by the decay due to density 
waves. Finally, as the disc parameters change with radius, 
the star eventually clears a gap in the disc and the dominant 
mode of its radial migration is switched once more. 

Figure |3 illustrates an interesting effect that emerges 
once the gravitational field of the disc is taken into ac- 
count. The disc introduces a non-spherical perturbation to 
the background gravitational field of the central mass. As 
a consequence, a Kozai-type mechanism causes occasional 
jumps of orbital parameters. A specific form of the disc grav- 
itational field is not particularly crucial for the onset of these 
hops, but the dissipative orbital decay and appropriate ar- 
rangement of time-scales are important (see Vokrouhlicky & 
Karas 1998 for more discussion). In fact, any non-spherical 
perturbation, for example general-relativity effects of frame- 
dragging near a rotating black hole, can potentially be re- 
sponsible for similar kind of oscillations. 

In the following sections we discuss integral properties 
of the whole cluster, and therefore we will sum over individ- 
ual stars. One should bear in mind that orbital resonances 
can set stars on very eccentric trajectory, increasing the cap- 
ture rate. Hence the presence of the disc may influence the 
overall cluster structure in its inner regions. 



2.4 The cluster 

Stationary state is maintained by a continuous inflow of 
fresh stars from the outer reservoir and this is accomplished 
by distinguishing between two populations of orbits: those 
which belong to the relaxed cluster with a prescribed dis- 
tribution fo and those which are influenced by the disc and 
belong to the inner cluster. The latter contribute by frac- 
tion /otinci/^r- The cluster as a whole is superposition of 
these two populations; the corresponding distribution will 
be denoted /* hereafter. Because the influence of the disc 
on stellar orbits is greater at small distances and likewise 
the effect of circularization increases inwards, therefore the 
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Figure 2. Long-term evolution of orbital elements for a single star. Different modes of star— disc interaction are distinguished by line 
type in the graphs, corresponding to the three subsequent phases of the mean motion (details in the text): in phase 1 the regime of direct 
collisions takes place (solid line); in phase 2 the regime of density waves occurs, in which a drops more rapidly (dashed line); finally, 
phase 3 describes migration of a star tidally coupled to the disc (dotted line). Time is scaled with respect to the initial dynamical period. 
To = 27r\/ GM, a^^^ ; semi-major axis is expressed in terms of the initial ao. 
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Figure 3. Similar to the previous figure, but for stars moving in superposed gravitational field of the central black hole mass and an 
axially symmetric thin disc. The spherical symmetry is now broken, bringing new features into the orbit evolution. Namely, sudden steps 
in the osculating elements occur at resonances between radial and vertical oscillation frequencies. Again, semi-major axis is scaled with 
respect to the initial value, ao, and time is given in units of To (only a part of the whole evolution is shown for better clarity). Short-term 
variations of the elements were filtered out in this graph. 



distinction between the two populations reflects also their 
spatial dimensions. 

First we introduce the distribution /o(ao, eo, xo, ojo), 
which characterizes the initial form of the cluster so far un- 
perturbed by the disc. It is determined by relaxation pro- 
cesses, acting on time-scale 



tr = 



G^ClnAM^n. 



(3) 



where a is the velocity dispersion, n, the number density 
of the stellar system, In A is Coulomb logarithm and C is a 
constant (in usual notation, C In A ~ 10^ ; e.g. Spitzer 1987) . 

If relaxation was neglected in comparison with the drag- 
ging of stars by the disc, one could write probability of flnd- 
ing a star on an orbit with semi-major axis a in the form 

/,(a; ao, eo, xo, cjo) da = dQ. /o(ao, eo,xo,iuo)tr^ ; (4) 

where J = ||da(a; ao, eo, Xo, (jJo)/dt|| and dfi is the phase 
space element. Meaning of this relation can be clarified by 
integrating eq. ||1J over a and obtaining the total number of 
stars in the inner cluster which have been injected from dQ. 
This number is equal to the fraction of stars in dQ multi- 
phed by the ratio of drag period to the relaxation period. 
At the same time, integration over the phase space of initial 
parameters gives the desired probability density in the inner 
cluster: 



/in(a)da= / dQo fo{ao,eo,xo,uJo)tr^J 



(5) 



The way in which the boundary zone is introduced between 
the inner and the outer populations is the main advantage 
(and limitation) of the adopted approach. It is assumed that 
SI is vacated solely by dissipative effects (interaction with the 
disc), while it is filled, at a constant rate, solely via relax- 
ation. In the numerical scheme this is achieved by requiring 
tinci to be shorter than relaxation time @ in the outer clus- 
ter, i.e. 



tin 



ao 



<t,. 



(6) 



where d is the rate of orbital decay (we substituted the value 
of d at initial time) . Another condition is imposed by requir- 
ing that the orbit intersects the disc, 



Rd > Ro 



ao (1 



eg) 



1 + eo cos ujo 



(7) 



A Monte Carlo code was used to evaluate the integral 
on the right-hand-side of eq. 0. We followed ~ 10^ or- 
bits with randomly generated starting values of osculating 
elements. Each time an orbit reaches a pre-determined value 
of tti (these are distributed in a logarithmical grid) the cor- 
responding bin of the probability-density array is increased 
by factor ld|~^. This contributes by 
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1 No sr^ \ da . , 

1 

to the distribution. Here, A^o is the total number of stars 
within the specified range of radii in the outer cluster, and 
A'c is the number of generated orbits intersecting the disc 
and satisfying conditions @ and Q . In typical calculation, 
A'c reaches a few percent of A'^. 

The code enables us to further distinguish two subsam- 
ples of the inner cluster. Stars on inclined orbits intersecting 
periodically the disc are counted to a dragged cluster. Once 
a star is embedded entirely into the plane of the disc and 
the mode of its radial migration is changed accordingly, we 
begin counting the star as a member of an embedded cluster. 
It turns out that these two subsamples populate different 
regions of the parameter space of the inner cluster, and so 
the distinction appears plausible. 

We have employed the BW distribution to populate the 
outer reservoir. Maxwellian velocity distribution is used as 
a boundary condition specifying the dispersion CTc at zero 
binding energy surface. As mentioned above, /o(ao) oc aj^* 
with ao being within (amin, dmax)- We set flmin = Rd and 
flmax ~ Rh, where -Rh is radius of the black hole gravitational 
dominance, 

7?h = GM. jal ~ 4.5M8cr2-oo pc, (9) 

written here in terms of (T200 = (Tc/(200 km s~^). Remaining 
orbital parameters, xa, eo and ljo, are distributed randomly 
with equal probability, ensuring spherical symmetry of the 
outer cluster. Corresponding spatial number density is 

n.(r)~(r/i?h no. (10) 

The velocity dispersion of the BW cluster scales roughly as 
a^^GMtjr. With help of eqs. @ and IIUL relaxation time 
JHJ can be rewritten as 

with ne = (no/lO'') pc^^. 

Next, we need to specify a relation between Uc and M, 
that would hold outside _Rh- We therefore adopt the empiri- 
cal relation (Ferrarese & Merritt 2000; Gebhardt et al. 2000), 
logj^Q Mg = a + /Slogj^o (7200 with a = 0.13 and /3 = 4 (cf. 
Tremaine et al. 2002). In our model, the log- log formula for 
crc(M,) is used to define the distribution of stars in the bulge 
overlaying our nuclear cluster and this way it unavoidably 
infiuences parameters derived from any observation. Other- 
wise, the forms of this relation in the two regions are quite 
independent of each other. Notice that simplicity is one of 
the reasons for adopting the BW outer cluster, but it is not 
necessary for the method as such. Likewise, different distri- 
butions of eccentricities can be accommodated. For exam- 
ple, we tested /*(e) oc e and /*(e) oc exp(e) without much 
impact on the results, i.e. on the final distribution of the 
system. 

The spatial size of the inner cluster is determined by the 
interplay between gravitational relaxation and dissipative 
effects in the disc. To estimate the size, two length-scales are 
naturally involved. Firstly, from eq. Q it is the disc radius 
i?d which limits the region of star-disc collisions. Secondly, 
we notice that t^(r) increases with radius at slower rate than 
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Figure 4. Typical examples of the structure of a cluster result- 
ing from its interaction with different models of an accretion disc. 
Top panels; distribution /»(a) of stars (solid line, in arbitrary 
normalisation) residing in a shell of given semi-axis (a in units of 
-Rg). Dashed, dotted, and dash-dotted lines stand for the three 
sub-components of the cluster, i.e. the embedded, dragged, and 
relaxed (outer) population. Bottom panels: mean orbital incli- 
nation (i) throughout the whole cluster, exhibiting the effect of 
flattening below lO^iig. 

iinci(a)- Hence, by substituting a— >r and setting t^ = find 
we define a characteristic radius 

For our exemplary disc models, one finds the following val- 
ues:^ (u) Ri„ = 4.73 x 10*7?g, (in) 4.36 x lO'^iig, (iv) 1.68 x 
10* -Rg, and (v) 9.80 x 10^J?g. Eq. (|T2l is useful for analyti- 
cal estimates of the cluster size. The smaller of the two radii 
J?i and _Rd corresponds to the size of the inner cluster. This 
estimate is in agreement with numerical simulations. 

Finally, we modified the above-described scheme and 
tested a model with the initial distribution /o being set to 
be Maxwellian in velocities. By comparing the results based 
on both approaches we found qualitative agreement between 
these two different formulations (see Karas et al. 2002 for 
details). 



3 RESULTS 

3.1 Radial structure of the modified cluster 

Let us first discuss the impact of star-disc interactions on 
the radial distribution of stars. Figure 2] provides connection 
between the properties of individual orbits and the overall 
evolution of the cluster. Parameters of the disc have been 
chosen identical as in Fig. |5| in particular, the outer edge 
7?d ~ 10^-Rg. Additional parameter is the stellar density, 

^ In the case (i) the estimate of drag time comes out more compli- 
cated, but we can safely consider that tr 2> tincl for any reasonable 
choice of model (i) parameters. 



Star-disc interactions in a galactic centre 7 



Disc (iii) 



Disc (ii)+(iv) 




2 3 4 5 6 1 2 3 
Logarithm of semi-major axis 

Figure 5. Ratio S{a\ fi) of two /« distributions obtained by con- 
sidering or neglecting the gravity of the disc in computations. Two 
cases are shown in each panel: First, the gravitational field of the 
disc was computed in consistent manner with the disc surface 
density SdC'") profile (solid line). Then, in order to demonstrate 
the dependency on the disc mass more clearly, we increased the 
11 parameter by factor of ten and recomputed the structure of 
the cluster with this new value (dashed line). The corresponding 
change of stellar density reaches order of several percent. 



no ~ lO^pc"''. In this case, radius Ri, defined by eq. 
is approximately 2 x lO^-Rg, i.e. twice the disc outer radius. 
Therefore, the transition from the inner to the outer cluster 
occurs at r = Rd- Clearly, the transition area introduced in 
this way is not a sharp value of radius - it depends espe- 
cially on eccentricity distribution of the stellar system, i?d 
of the disc and tinci/ir of the both components. The resul- 
tant number density, /, (a), is closely matched by a broken 
power-law with a slope depending on the dominant regime 
of star-disc interaction. Notice that the asymptotic form is 
/, oc a^^* for large a. 

The impact of disc gravity can be demonstrated by 
"switching" the gravitational field of the disc on and off 
in computations. This is achieved by setting the value of 
H parameter and computing the corresponding f,{a;jj,). In 
Figure 1^1 we show resulting ratios 



S{a;fi) 



/. (a; fi^O) 



(13) 



of two such distributions. Each curve was obtained with 
two identical models, except for the inclusion or omission of 
the disc gravity. The two frames differ only by the adopted 
model of the disc. It turns out that perturbation to the stel- 
lar orbits by non-zero disc mass speeds up the orbital decay 
and, consequently, it decreases the number density of stars 
in the dragged cluster. This is visible even if the mass ratio 
is small, i.e. of the order of ^ ~ 10~^ 10~^. Naturally, 
the influence of the disc gravity increases with jj, increasing. 
Details of profile of the gravitational field also matter: we 
see, for example, that the disc (iii) affects /, more than the 
disc (ii)-l-(iv) if their mass parameter /i is comparable. The 
reason of this behaviour is that the former case has a fiat- 
ter curve of Sd('') (cp. Fig. Q extending farther out from 
the black hole than the latter case, and hence it exhibits 
stronger impact on the cluster distribution. 

Furthermore, we have included the orbital decay by 
gravitational radiation. Indeed, it is highly desirable to 
estimate this type of effects also because of forthcoming 
gravitational-wave experiments. We employed a simple ap- 
proximation (Peters & Mathews 1963) and we found that 
gravitational radiation infiuences the final /, (a) only very 



near the centre (see Subr & Karas 1999 for details) . For the 
discs considered here, this occurs typically below < 10-Rg 
and, therefore, it does not leave any remarkable imprint on 
the results. Notice that this conclusion does not contradict 
the results of Narayan (2000) who considered case of low- 
density fiows, and therefore he found the gravitational radi- 
ation to be more important for the orbital decay than the 
hydrodynamical dissipation, also at larger distance from the 
centre. 

Because the system is non-spherical but axially sym- 
metric and stationary, it is natural to portray its structure in 
a two-dimensional spatial projection. Figure |S| shows differ- 
ent sections of the cluster arranged in four columns. Starting 
from the left side, columns (a) and (b) represent the spatial 
number density n, (r) in the meridional plane across the clus- 
ter. Individual frames are centered on the black hole. Fur- 
ther, two columns on the right show the projected (column) 
density, i.e. the total number of stars integrated along the 
observer line of sight. Two different lines of sight were cho- 
sen, perpendicular to (c) and co-linear with (d) the symme- 
try axis. Interesting morphological patterns can be revealed 
by comparing different views. 

Let us describe the meaning of individual rows of Fig.|S| 
In the first row, the density profile of the unperturbed BW 
distribution is plotted. It corresponds to the outer cluster 
which is clearly spherical and serves as a starting config- 
uration for our computations.^ We show this case for the 
sake of comparison with the structure of modified clusters 
in subsequent rows. Again, we use the disc models (iii) and 
(ii)-l-(iv) as two representative cases. 

Comparing different representations of the cluster one 
can clearly observe the impact that star-disc collisions have 
on the cluster structure. In particular, it is the increasing 
oblateness of the stellar population in the core and, in some 
cases, the tendency to form an annulus of stars. The effect 
is best seen in the bottom row in Fig. |S| However, details 
of the modified cluster structure obviously depend on the 
type of accretion disc that we choose. The reason of dif- 
ferent structures stems from the different modes of radial 
transport of the stars from the cluster to the central hole. 
In the following paragraphs we examine several symptoms 
of the growing deformation of the cluster. 



3.2 Anisotropy 

The rate of dragging the stars by the disc is higher for 
counter-rotating stellar orbits than for co-rotating or aligned 
ones, hence, it is natural that the initial state of isotropy 
and sphericity becomes gradually violated. Less obvious is 
the fact that toroidal structures may be formed as a conse- 
quence. To see the effect we inspect graphs of mean incli- 
nation {i)(a) (bottom row of Fig. 0. The above-mentioned 
populations of cluster members are clearly distinguished, in 



^ Notice that n,(r) exceeds 1.7 X 10^^ pc~'^ below O.liJd. It is 
thus out of density scale in the detailed plot of the central region 
- column (b), top frame. However, such a high concentration is 
only a formal artefact of the BW model in which stellar density 
diverges in the origin, n,{r) oc r"^/^. The divergence does not 
occur in our cluster interacting with the disc, in which case the 
BW solution is only used to describe the outer reservoir. 
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Figure 6. The spatial density n,(r) and the corresponding projected density of the cluster arc shown using logarithmically spaced 
levels of shading. Columns (a) and (b) represent the meridional section at two different scales, namely, 4iJ(j and 0.2_R;j across (radii are 
expressed in terms of the disc outer radius, Rd)- Next columns are the edge-on (c) and the face-on (d) projections of the cluster. Across 
columns, the upper row shows the referential cluster (n* oc r~'^/'^). In subsequent rows, the system has been already modified via the 
interaction with two types of discs, case (iii) and case (ii)-l-(iv), as indicated. 



particular, small values of (i)(ffl) indicate a flattened struc- 
ture of the cluster with a population of disc co-rotating 
orbits at small a. While the outer cluster approaches an 
isotropic system ((i) 7r/2) for large a, the dragged cluster 
is clearly flattened to roughly a constant value of mean incli- 
nation, ((i) ~ 0.2). It drops further in the embedded cluster 
where {i) — > 0. This behaviour provides another substanti- 
ation for distinguishing different populations that form the 
cluster. 

The limiting value of moan inclination can be derived by 
the following simple analytical argument (which agrees with 
the results of computations). Probability density of finding 
a body in (a, a + da) x {x,x + Ax) is 



f*{a,x) = / dao da;o di/o dwo 



fo{ao,xo,yo,oJo) 
tr(oo) 



da 
dt 



X 6{x-x{a;ao,xo,yo,i^o)), 



(14) 



where integration is performed over the whole phase space 
of initial values. The mean inclination is 



(.)(a)./-(a)/d.arccos./,(a,.), 
where the normalization factor /« (a) is 

f,(a) = / da;/*(a,a;). 



(15) 



(16) 



The mean inclination can be estimated analytically by 
neglecting gravity of the disc and assuming all orbits to be 
initially circular with random inclination, 



/o(ao,xo,yo,wo) = /Bw(ao) (5(j/o - 1) ■ 



(17) 
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Figure 7. Top; eccentricity distribution /« (e) in the inner cluster. 
Bottom: ratio S{e; fi) showing the effect of disc gravity according 
to eq. 1211 . Parameters of the model are described in the text. 



With the aid of eq. the integral in eq. H15^ can be esti- 
mated in the form 



{i)(a) 



2\/2 7r 



21n8-21n(l-a:;n,ax)-3 + 3;n 



(18) 



(a^max denotes cosine of limiting inclination that distin- 
guishes stars of the dragged and embedded cluster). For thin 
disc models, 1 — Xmax ~ ^(H/r)^ is roughly constant, reach- 
ing values of the order of 10"^ -i- 10^^. Hence, we obtain 



{i){a) 



V2n 



ln(l — Xn 



' 0.3, 



(19) 



which only slightly overestimates values we found by Monte 
Carlo computations. 

In the velocity space, anisotropy can be measured by 
means of factor 



/3a 



2{v?} ■ 



(20) 



which is defined in terms of transverse and radial velocity 
components, v-r and v^, averaged over the cluster. Result- 
ing values of |/3a| are typically of the order of unity, zero 
value corresponding to isotropic distribution. Averaged over 
the whole cluster, our models give /3a ~ —1.5, while values 
within interval (—5, —4) are obtained for the inner cluster. 



3.3 Diversification of eccentricities 

By gradual circularization, every captured star contributes 
to the resulting profile /« (e) de proportionally to time it 
spends within (e, e+de). The orbital evolution slows-down at 
the late phase and this leads to higher abundance of nearly 
circular orbits. Prevailing circularisation is therefore antic- 
ipated because of gradual energy losses and the resulting 
orbital decay. Fraction of highly eccentric orbits in the clus- 
ter is small, yet it is important to know because these orbits 
bring stars close to the black hole where they can be prefer- 
entially captured or destroyed. The effect of non-sphericity 



of the gravitational field complicates the long-term tendency 
because it can lead to occasional increases of eccentricity. 

We can distinguish two concurrent mechanisms that de- 
termine eccentricities. Firstly, the selection rules JSJ and 
define a loss cone in the space of initial osculating elements. 
Stars on eccentric orbits are continuously shipped to the 
inner cluster and this mechanism accentuates the propor- 
tion of highly eccentric orbits with pericentres below _Rd, 
hence, such orbits are more abundant in the final distribu- 
tion. Secondly, it is the interaction of captured stars with the 
disc which shapes the inner distribution. Combining the two 
mechanisms, we obtain a double-peaked distribution with lo- 
cal maxima at low and at high eccentricities and depression 
in the middle (see Figure |7|l. This structure corresponds to 
an intermediate phase of the cluster evolution which is seen 
in simulations of Ranch (1995; cp. his Fig. 8). 

When the disc gravity is taken into account during the 
orbit integration, resonances lead to eccentricity jumps, this 
way contributing to highly eccentric orbits. We examined 
the overall impact on the distribution function by means of 
ratio 



S{e-fj.) 



(21) 



/*(e;M = 0) 

analogically to eq. 11311 . Figure |7| suggests that a non- 
spherical perturbation to the gravitational field leads to 
higher abundance of more eccentric orbits. 

3.4 The mass stratification 

Let us now examine the mass spectrum of the cluster. We as- 
sume a "canonical" form of the initial mass function J^o(A'It) 
(Salpeter 1955) together with the empirical relation for 
E,(M,) of main sequence stars (e.g. Binney & Merrifield 
1998; Lang 1980): 



AAo(M,) oc M- 



s*(M*) oc m; 



(22) 



Stellar masses are generated in the range 0.5Mq < M, < 
IOMq. Notice that the assumption 12211 is needed for defi- 
niteness of the following example, but its special form and 
the assumed range of masses are not required by the adopted 
approach. This way we obtain profiles resembling, typically, 
those for a cluster which is continuously supplied by equal- 
mass stars at the lower boundary of the mass spectrum. The 
resulting distribution follows partly from the fact that light 
stars represent a major fraction of inserted bodies, partly it 
is due to their relatively slow migration towards the centre. 
Also, radial transport of light bodies in the disc is less effi- 
cient than that of heavy ones, except for very massive stars 
which tend to open a gap in the disc. Once this happens, 
a qualitatively different mode of radial motion takes place, 
typically orders of magnitude slower than under the regime 
of density waves excitation. 

Gradual change of N{Mt) arises here due to unequal 
rate of orbital decay of stars with different masses. The effect 
of mass stratification can be estimated by 



AA(M.) oc AAo(M.) \h{M,) 



(23) 



See Figure|H|for profiles resulting from computations (values 
of the best-fit power-law index are given in the plot). Dif- 
ferent curves show the stellar-mass distribution within the 
corresponding radial distance from the centre. 
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Figure 8. An exemplary graph of the mass function (with arbi- 
trary normalization) for a cluster modified by the interaction of 
its members with diff'erent types of discs: case (iii) is shown on 
the left, case (ii)+(iv) on the right. Individual curves represent 
the mass function of a sample of stars occurring below a given ra- 
dius: lO^Rg (solid), lO^Rg (dashed), lO^R^ (dotted) and lO^Rg 
(dash-dotted) . Numbers are given with the curves to indicate the 
corresponding power-law indices of best-fit lines. 
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Figure 9. Velocity profile along the line of sight of the in- 
ner cluster, integrated across a column of cross-sectional radius 
Rs = 10-*iJg (left panel) and = lO^Rg (right panel). SoUd 
and dashed lines represent different inclinations of the observer: 
5 = 0° and 60° , respectively. Growing anisotropy of the modified 
cluster produces the dependence of measured line profile on ob- 
server's view angle, i.e. 1=1 (vz \ §). The thin dotted line stands for 
the referential BW distribution, which is spherically symmetric. 



Let us discuss case of the disc (ii)-l-(iv) in Fig.|Hlin more 
detail. Majority of stars in the region below 10^i?g belong 
to the outer cluster. The mass function remains almost un- 
modified in this region. However, selecting only stars inside 
10* Tig into consideration, substantial fraction of orbits are 
aheady aligned with the disc. Their radial velocity of migra- 
tion is proportional to EJ^ and, through this dependence, 
it is cx Mf'^. Hence, employing eq. 1231 we estimate that 
the mass function index decreases by 0.6 and it reaches the 
value —2.95. Going still closer to the centre, we find more 
stars that are fully embedded in the disc and migrate in 
the regime of density waves, i.e. with radial velocity oc M, . 
They cause a further drop of the mass index and the ex- 
pected value in this case is —3.35. This trend is in good 
agreement with the behaviour we see in numerical simula- 
tions. 

A modification of the power-law profile occurs at the 
upper boundary of the stellar mass range as a consequence 
of changing the mode of radial migration. This is evident 
in the case (iii). Considering r<10"^i?g, majority of stars 
again belong to the embedded cluster, but massive stars 
now succeed to open a gap in the disc and, therefore, they 
continue migrating to the centre with the radial velocities 
independent of M* . For this reason the mass-function index 
is raised closer to its initial value. Within the innermost 
region (below 100 iZg) the limiting mass for the gap opening 
decreases and the transition between two different slopes 
moves to lower M, . 

3.5 Velocity dispersion 

Let us examine integral properties of the cluster that refiect 
the state of the inner region in terms of potentially observ- 
able quantities: the shape of a synthetic spectral line I{vz) 
(i.e. intensity as a function of line-of-sight velocity Vz near 
the projected center of the cluster) and velocity dispersion 
(Jz along line of sight. It is worth noting that similar spectral 
profiles were computed for the reference cluster already in 
the original Bahcall & Wolf (1976) paper, so one can check 
how the picture is changed by the presence of the assumed 
disc and what imprints could be traced in the integrated 
lines from the nucleus unresolved to individual stars. 



In Figure El the left panel shows the line-of-sight pro- 
file that was obtained by averaging over a column of radius 
Ra = 10*i?g, centered on the cluster core. This profile de- 
termines an infinitesimally narrow line (in a local frame of 
each star). We neglected all other effects that might influ- 
ence the final spectral feature and assumed the disc (iii) with 
Ri = 10''7?g. Another view of the same system is shown in 
the right panel, where we chose lower spatial resolution for 
comparison. Also in these profiles one can recognize the ten- 
dency of the stellar sub-system to be fiattened. Stellar ver- 
tical velocity component decreases with respect to the com- 
ponent parallel with the equatorial plane, which manifests 
itself by gradual decrease of the line width as the view angle 
becomes parallel with the symmetry axis. Local maximum 
of the line occurs around v ~ sin(^)iiK(-Rd)- In some cases, 
this secondary peak may exceed the central maximum and 
dominate the whole profile. High-velocity tails of the line 
profiles are also noticeably affected. 

In Figure lTUl we show predicted dispersion using a cir- 
cular aperture of radius R^. In the case of initial BW distri- 

1/2 

bution, az scales roughly proportionally to M, . The slope 
is somewhat modified when the range of semi-major axes is 
restricted in the reference cluster (in particular, we assumed 
Rd ^0,0^ Rh)- If the cluster is modified by interaction with 
the disc, a threshold mass can be found in the az{M,) depen- 
dence. Its value depends on the fraction of stars belonging to 
the inner cluster. At the lower end of the M, range, disper- 
sion resembles the previous unperturbed case, because the 
anisotropic inner cluster is entirely concealed by the outer 
one. On the other side, for large M, we observe a steeper 
slope of the cJz{M,) relation - a consequence of the newly 
established radial density profile n«(r). The exact form of 
this relation depends on the view angle, which is another 
manifestation of the inner cluster anisotropy. 

3.6 Mass inflow from the outer cluster 

The rate of mass infiow in stars, Ms, characterizes their 
radial drift from the outer reservoir. Related to Mg is the 
induced energy dissipation Kb, also caused by star-disc col- 
lisions. These two quantities can be expressed in a common 
way. To this aim we introduce a new variable, A, denoting a 



Star-disc interactions in a galactic centre 11 




2.5 



rate M of the gas in the disc. As for Ks, we expressed this 
quantity as difference Vi — Vi of the initial and final potential 
energy of a star, where V{ was estimated from eq. Q. 

We remind the reader that a uniform distribution of 
initial eccentricities was assumed, but this limitation is not 
very important because only a narrow range of large initial 
eccentricities populate the resulting distribution.^ Eqs. i27l 
and ^ can therefore serve as useful order-of-magnitude 
estimates. 
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Figure 10. Velocity dispersion CTz along the line of sight as a 
function of the central mass M, (top panels). Different line- types 
represent observer's view angle: § = 0° (solid), 30° (dashed) and 
90° (dash-dotted). For the sake of comparison, thin dotted line 
is an unperturbed profile of the outer cluster. In bottom panels 
we magnify the differences showing the ratio of Cz with respect 
to the BW case. 



quantity that is formed by contributions carried by individ- 
ual stars, as they sink to the centre {Ks or Ms are just two 
examples of such quantity). The corresponding flow rate is 

/•Omax /•! /"l p'^ 

A = dao / deo / da::o / da;o .4 /Bw(ao, eo) (ao) 

•^a„iin Jo J -1 Jo 

X e(-ao - ao/tr) e(i?d - Ro), (24) 

where the two theta-functions limit the region of parameter 
space from which the stars are captured by the disc; cf. 
conditions Q and iQ. We obtain Ms by setting 

(25) 

in the integrand of eq. (I24II . Likewise we obtain Ks for 

A-^ AE = GM.M, (— - —) . (26) 
V flf aoJ 

In both cases we performed integration numerically. One can 
omit the leading theta-function in eq. 12411 because the drag 
time is shorter than relaxation time for majority of initial 
conditions (we assumed customary values of accretion disc 
parameters, Mo.i = 1, ao.i = 1). Even if relaxation cannot 
be immediately neglected, we may still ignore the mentioned 
term in eq. II24II and substitute R^^Ri. 

Further, we write analytical estimates for Ms and Ks, 



Ms = 1.7 X IQ-^ Ml^'^nl 



ks^lQ-^Ml^^nl 



Mq 





' Ri 




WR 


3 

£ 


' MqC 





(27) 



(28) 



where £ and J- are order-of-unity and slowly varying func- 
tions for 0.01 7?d < flmin ^ -Rd, Rd flmax < lOOiid . 

The total accretion rate onto the black hole is a sum of 
Ms, which involves stars from the cluster, and the accretion 



4 DISCUSSION AND CONCLUSIONS 

We have developed the model of a cluster departing from 
its initial sphericity and isotropy as a result of star-disc 
collisions. One can anticipate various observational conse- 
quences of this interaction, and several of them were ex- 
amined in Sec. |3 Hereafter we list other implications but 
we defer detailed comparisons with actual evidence until a 
more complete picture is available. Let us remark that ob- 
servations can be confronted in two areas: firstly, it is the 
overall form of the stellar distribution in the nuclear clus- 
ter as discussed above (which however corresponds to sub- 
parsec scales, not yet sufficiently well resolved), and then it 
is a discussion of individual events of star-disc collisions or 
disruptions, whose frequency is linked with the rate of or- 
bital decay and which may be detected by a sudden releases 
of X-rays or by gravitational-waves emission. 

The orbital decay of stars near a black hole is indeed 
relevant for forthcoming gravitational wave experiments, be- 
cause the gas-dynamical drag needs to be taken into account 
in order to predict waveforms with sufficient accuracy. It has 
been estimated (Narayan 2000; Glampedakis & Kennefick 
2002) that this effect can be safely ignored at late stages, 
shortly before the star plunges into the hole, if accretion 
takes place in the mode of a very diluted flow. However, the 
situation is quite different in case of AGN hosting rather 
dense nuclear discs (Subr & Karas 1999). Consistent eval- 
uation of gravitational radiation from dense stellar systems 
would require reformulation of our model within the gen- 
eral relativity framework, which is beyond the scope of this 
paper, but we estimate that the motion of only a minor 
fraction of stars very near the centre should be noticeably 
affected by the emission of gravitational waves. One can ex- 
pect that exact calculation would lead to even a weaker ef- 
fect because our stellar system tends to a flattened rotating 
configuration that does not radiate in the limit of perfect 
axial symmetry (luminosity emitted in gravitational radia- 
tion by particles orbiting in a ring decreases exponentially 
with the number of particles; Nakamura & Oohara 1983). 
In other words, gas-dynamical effects most likely dominate 
over the gravitational radiation as far as their influence on 
the cluster structure is concerned. 

The modifled cluster structure is obviously pertinent 
for various studies of the black-hole feeding problem and, 



^ For another form of initial eccentricity distribution, an esti- 
mate of the fiow of mass and kinetic energy can be obtained by 
multiplying formulae 1271 and I28i by factor /*(eo = 1). We 
verified accuracy of this scaling in two cases, /*(eo) = 2eo and 
/, (eo) = 0.58exp(eo), and we found that it indeed agrees with 
numerical computations. 
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vice- versa, for the issue of a possible feedback that a super- 
massive black hole may exhibit on the the host galaxy bulge. 
However, the well-known empirical relations between the 
black hole mass and the bulge properties apply to the region 
substantially exceeding Typically, M,/Mbuigc ~ 10~^, 
while the two masses should be comparable where the clus- 
ter structure is to be directly influenced by the central hole, 
which is what we assumed here. An interesting direction in 
which our model can be further developed therefore con- 
cerns a more realistic form of the outer cluster and the way 
it is connected with the galaxy bulge. This could help to 
build a model for the ac{M,) relation extending beyond the 
domain of black hole gravitational influence. In this con- 
text it is worth to mention a recent idea (Miralda-Escude & 
Kollmeier 2004) that star-disc collisions can indeed be seen 
as a self-regulating process that helps to feed the central 
black hole and control its growth. 

Inherent limitations persist in our present discussion, 
namely, we have not incorporated a fully self-consistent 
treatment of the disc gravity. Even though we computed 
the gravitational fleld across a sufficiently large domain of 
space, we did not account for the feedback which star-disc 
collisions exert on the disc structure. Introducing some kind 
of a clumpy model of the disc will be very interesting, as 
it may exert a more substantial effect on the cluster struc- 
ture by elevating the impact of star-disc collisions. Different 
configurations of a clumpy disc have been proposed (Kumar 
1999; Fukuda, Habe & Wada 2000; HartnoU & Blackman 
2001; Kitabatake & Fukue 2003) which will be relevant for 
this investigation. 

Also relevant is the idea of enhanced star formation 
caused by supernovae that are triggered in a self-gravitating 
disc, changing substantially its structure (Collin & Zahn 
1999). We can expect that the inner cluster will be af- 
fected in different way in those galactic nuclei where nuclear 
starbursts occur simultaneously with AGN phenomenon 
(Tenorio-Tagie et al. 2003; Gonzales Delgado et al. 2004). 
The circumnuclear starburst phenomenon has been indeed 
indicated in some AGN on ~ lO'^ pc scales, which could cor- 
respond to the outskirts of a dusty torus surrounding the 
central black hole (e.g. Aretxaga et al. 2001; Heckmann et 
al. 1997; Schinnerer et al. 2000). Star formation could then 
take place at outer parts of the nuclear cluster. 
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